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Abstract 

The change in orientation in generation space (rotation) of the 
fermion mass matrix with changing scales can lead to flavour- violations 
through just the kinematics of a non-diagonal mass matrix. Such ef- 
fects for the reactions: e^e~ — > e^/i^, e^r^, //^r^, and for the de- 
cays of vector bosons into the same channels, are calculated following 
a method suggested earlier which gives the differential cross section 
for each reaction and the branching ratio for each decay mode in terms 
of an overall normalization depending only on the speed at which the 
mass matrix rotates. A rotation speed estimated earlier, under cer- 
tain assumptions from the fermion mixing angles and mass ratios, is 
found to give the above effects at a level readily detectable in mod- 
ern high sensitivity experiments such as Bepc, Cleo, BaBar and Belle, 
at least in principle. The observation of these effects would not only 
confirm the concept of a rotating mass matrix with a significance on 
par with the running coupling constant, but also offer some valuable 
insight into the origin of fermion generations. However, a negative 
result cannot unfortunately rule out the rotating mass matrix since 
the effects deduced here from this mechanism alone could in principle 
be cancelled by other rotation effects. 



1 Introduction 



With the running couphng constant now a famihar concept aheady amply 
verified in experiment, one would not be surprised that the fermion mass 
matrix also varies with changing energy scales. That its eigenvalues, namely 
the fermion masses, actually do run has already been experimentally ver- 
ified in certain circumstances What we consider in this paper is the 
scenario when the fermion mass matrix changes also its orientation (rotates) 
in generation space as the scale changes. 

Theoretically, of course, there are good reasons to expect that the fermion 
mass matrix will rotate with changing scales. Even in the Standard Model 
mass matrix rotation occurs so long as there is nontrivial mixing between 
the up and down fermion states. The renormalization group equation there 
satisfied by the mass matrix U of the up fermions pf: 

contains already at leading order a term on the right which is nondiagonal in 
the eigenstates of U when the mass matrix D of the down fermions is related 
to f/ by a nontrivial mixing matrix, so that a U matrix diagonalized at one 
scale can no longer remain diagonal at another scale, or in other words, it 
will rotate as claimed. The same argument holds for the mass matrix D for 
the down fermions. This fact was pointed out for quarks already long ago 
in, for example, [Q, and now that nontrivial mixing for leptons has been 
strongly indicated |^, |, ^, H, if not already confirmed by experiment, 
the same conclusion can be drawn for leptons also. Next, looking beyond 
the present Standard Model, one encounters further possible mechanisms for 
driving the mass matrix rotation. Indeed, the mere fact that different gener- 
ations of fermions can rotate into one another, as the mass matrix rotation 
implies, already means that they are not distinct entities as once conceived 
but just different manifestations of the same object, like the different colours 



of a quark, related presumably by some "horizontal" symmetry [T^. This 
suggests new forces which can change the generation index, and hence drive 
directly mass matrix rotations, in addition and in contrast to the "indirect" 
mechanism in (|1|) via mixing. 

However, even apart from theoretical prejudices as above indicated, there 
are, in our opinion, already some empirical indications which argue strongly, 
though perhaps as yet only circumstantially, in favour of fermion mass matrix 
rotation. These come about as follows. As is well-known, quarks and leptons 
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exhibit remarkable mass and mixing patterns. First, both their mass spectra 
are hierarchical, with the masses falling by one to two orders in magnitude 
from generation to generation. Secondly, the mixing matrices (i.e. CKM 
TT|, O] for quarks and MNS O] for leptons) which parametrize the relative 



orientations between up and down states seem roughly similar in shape for 
quarks and leptons, at least as far as is already known about the latter 
from the data in neutrino oscillation 0, |], |^, 0, |, only with the off- 
diagonal elements generally much larger for leptons than for quarks . Neither 
of these features have any explanation in the present Standard Model, in 
which they are just taken for granted, and between them they account for 
some two-thirds of the Standard Model's 20-odd empirical ("fundamental") 
parameters. However, as was pointed out in a recent note |jl5[, if one assumes 
that the fermion mass matrix rotates at a certain speed, then these features 
can all be very simply understood. 

That this is the case can be summarized as follows. Once the mass matrix 
is allowed to rotate with changing scales, then the usual definition of fermion 
flavour states as its eigenstates will have to be refined since, the eigenstates 
being now scale-dependent, one has to specify at what scale (s) the fermion 
flavour states are taken as the eigenstates. If the fermion masses are defined 
as the eigenvalues of the mass matrix at the scales equal in value to the 
masses themselves, as is usually done, then it would be natural to define the 
corresponding fermion flavour states as the eigenstates at the same scales 
also. It follows then that the state vectors of, say, the t and b quarks will be 
defined at different scales, so that even if the U- and D-quark mass matrices 
share always the same orientation at the same scale, the rotation of the mass 
matrix from the scale of the t mass to that of the b mass will already imply 
a disorientation between the t and b state vectors, or in other words a CKM 
matrix element Vtb represented by the direction cosine between the t and b 
vectors which is different from unity. More generally, it can be seen along 
the same lines that a rotating mass matrix will generate not only a nontrivial 
mixing (CKM or MNS) matrix but also nonzero masses for lower generation 
fermions even when one starts with neither. Indeed, what was shown in 



was that all existing empirical information on the fermion mass hierarchy 
and mixing pattern, (excepting for the moment only CP- violation) can now 
be understood as consequences a rotating mass matrix in the above manner, 
at least qualitatively but in some cases even quantitatively. We regard this 
as a rather strong though indirect empirical indication in favour of fermion 
mass matrix rotation. Indeed, if this interpretation of fermion mixing and 
mass hierarchy is accepted, then it implies a rotation speed considerably 
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greater than that driven by the mechanism (|T|) in the current Standard Model 
framework and suggests a driving mechanism from beyond that. 

Given these indications, both empirical and theoretical, it would be nat- 
ural to enquire what other physical implications a rotating fermion mass 
matrix may have which can be tested directly by experiment. One obvious 
candidate is flavour-violation, for a rotating mass matrix will not remain di- 
agonal in the flavour states at scales other than the scale (s) at which these 
flavour states are defined. And since reaction amplitudes depend in general 
on the fermion mass matrices, these too can become flavour-nondiagonal 
leading thus to flavour-violating reactions. This possibility has already been 
considered in general terms in, for example, [|1^ in some detail. We need 
here only to give an outline of its physical significance. 

The importance of flavour-conservation as a possible fundamental con- 
cept, of course, has long been recognized and its consequences subjected to 
rigorous experimental tests. In particular, the very stringent bounds set on 
/i — > 67 and fj, — > eee decays, which are at present respectively 1.2 x 10~^^ 
and 1.0 X 10~^^ for the branching ratio over the predominant, but already 
weak, decay for /i, tend to give the impression that any violation of lep- 
ton flavour would have to be extremely small. However, such a conclusion 
may be premature for not having taken account of the possibility that the 
lepton mass matrix rotates with the energy scale. If the mass matrix does 
rotate, then the fact that flavour-conservation has been stringently tested in 
^ — > 67 and fi — ^ eee at the fi mass scale, where the state is by defini- 
tion diagonal, is by itself no guarantee that flavour- violation will be equally 
small in another reaction at another scale where the mass matrix may have 
rotated to another orientation so that the /i state is no longer diagonal there. 
Such flavour-violating effects can in principle occur by virtue of the mass ma- 
trix rotation even when there are no explicit flavour-changing neutral current 
(FCNC) couplings in the action, and can be sizeable even when FCNC effects 
are small. We have therefore suggested for them the term "transmutation" 
in for distinction, which we shall adopt also in this paper. The size of 



transmutation effects, if any are observed, and their variation with energy 
will give indications on how the rotation of the mass matrix is driven, the 
knowledge of which may in turn shed light on the origin of generations itself, 
a basic question in particle physics that has already been with us for many 
years. For this reason, we suggest that flavour-violation be routinely tested 
in experiment whenever conditions are favourable. 

For pursuing this program, a method was developed in |r7| for calculating 
the flavour- violation effects due just to the kinematics of a rotating mass ma- 
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trix. In this, one treats the rotating mass matrix as a proposition in isolation 
without enquiring from what mechanism this rotation originates and with- 
out taking account of other possible rotation effects which might in principle 
accompany the mass matrix rotation. Such a preliminary attitude, we think, 
is reasonable given that the only empirical evidence one has so far is for the 
rotating mass matrix alone with no direct hints yet of the mechanism 
driving it. The result of such a calculation will serve to gauge what size 
flavour-violation effects might in principle be expected from a mass matrix 
rotating at a given speed with the aim of providing an indicator for experi- 
menters planning an analysis along these line. However, it is not to be taken 
as a necessary prediction of a rotating mass matrix under general circum- 
stances for the following reason. Reaction amplitudes depend on quantities 
other than the fermion mass matrix, such as, say, interaction vertices, which 
may in principle also rotate, and any theoretical mechanism for driving the 
mass matrix rotation may also imply vertex rotations. And these other rota- 
tional effects may modify or even cancel the effects computed by the above 
method from the kinematics of the rotating mass matrix alone. Indeed, in a 
detailed calculation performed specifically in the Dualized Standard Model 
(DSM) framework that we ourselves suggeted which is reported in a separate 
paper |]T^, such a cancellation is found in fact to occur. Hence, calculations 



done with the method of |T6|, [T^ as those in the the present paper have to 
be taken with this reservation in mind, which reservation was unfortunately 
not made clear in the earlier references because it was not clear then even to 
ourselves. 

In the present paper, we choose to investigate transmutation effects in the 
following reactions by the method suggested in |]l^| which, though originally 



developed for photo-transmutation, can be adapted to the present case with 
but minor modifications: 

e+e" — > e+r", r+e"; (2) 

e+e" — i> e+/i", /i+e"; (3) 

e+e~ — > jJi^T^ .T^jjT. (4) 

The obvious practical reason for investigating these reactions is that there are 
several high intensity machines in operation, such as BEPC, CESR (Cleo), 
PEP II (BaBar), and KEK II (Belle), which appear capable of observing 
these flavour-violating effects to high accuracy, besides LEP, which though 
now turned off, has left still masses of data which could be useful for the 
same purpose. 
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Lepton transmutation in reactions (j^) — can proceed by, for example, 
the processes represented by the Feynman diagrams in Figure |T[ At the 
energy at which an experiment is performed, the amphtudes for these pro- 
cesses, being dependent on the lepton masses, are diagonal in the eigenstates 
j = 1, 2, 3 of the mass matrix at that scale but not in general, by the reason- 
ing above, diagonal in the flavour states e, /i, r. And this fact alone could be 
enough to give lepton transmutation as a result. The transmutation reactions 
0) — ® can of course occur also via other processes, such as higher order 
photon-exchange diagrams or Zg-exchange, but the effects from these for the 
energy range of present interest are small and will be neglected. What can 
give sizeable contributions, however, is the formation of vector bosons in the 
intermediate state followed by their subsequent (transmutational) decays, as 
represented by Figure || (b). 




Figure 1: Feynman diagrams for the transmutation amplitude. 
Figures ^ and |^, then, are the only processes we shall consider in this pa- 



per, where we shall show that, by following the procedure suggested in ||T7 
one can calculate explicitly the differential cross sections for the 3 transmu- 
tational reactions (H), @, and (^), given any rotating mass matrix. The 
rotating mass matrix itself will figure only in the normalization of the cross 
sections, not in their angular or spin dependence both of which are given 
essentially just by kinematics. 

As a numerical example for the sort of cross sections one might expect for 
these transmutation reactions, let us consider in particular the reaction (^. 
The normalization of the cross section of this reaction is given by the rotation 
angle between the /i and r states at the energy scale of the experiment, a 
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(a) 



Figure 2: Decay amplitudes 



good estimate for which can be read aheady from the existing data on fermion 



mass and mixing patterns interpreted as rotation effects as in |15], say from 



Figure 3 of that paper. Specifically, for ^/s = 10.58 GeV at the mass of the 
T{AS) at which BaBar is run, one obtains in this way an estimate for the 
integrated cross section of reaction (|^) of around 80 fb which translates to 
as many as a thousand events in the data they have already collected from 
their first year's run of 20 fb~^, assuming 100 percent detection efficiency, 
and if so should be readily detectable. 

However, as already stressed above, these estimates for transmutation 
cross sections are not to be interpreted as definitive predictions for the quoted 
rotating mass matrix since they can be modified by other rotation effects. 
Nevertheless, they are of interest in giving an idea of the size of flavour- 
violation effects that a rotating mass matrix can in principle generate. Given 
the smallness of flavour- violation at the fi mass scale in /i — >■ 67 and fi eee 
decay, any detection of flavour- violation in the reactions (0) — (^) at BaBar 
or similar experiment at a different scale would be a positive indication for 
mass matrix rotation, although a negative result, by virtue of the preceding 
observation, would not be able at present to rule it out. 



2 The Reaction Amplitude (a) 

Consider flrst the one-photon exchange diagram of Figure |l](a), which will 
be seen to give the dominant contribution to the two reactions (0) and @. 
It does not contribute to the reaction (^), which can proceed by one-photon 
exchange in e^e~ collision only when both and e~ transmute, but this will 
be so far down in magnitude as to be negligible for present consideration. 
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Then according to the procedure suggested in at any given energy ^/s, 
the transmutation amphtude for the reaction : 

e+C — (5) 
is given by a rotation in generation space, thus: 

M^'^^ = Y,SlMfs^,, (6) 
j 

from the diagonal amphtudes ■M.'^f^ for the reaction: 

e+£- e+ij (7) 

for the mass eigenstate j at the scale ^/s with eigenvalue mj, where Saj = 
{j\a) is the rotation matrix in generation space which relates the triad of 
lepton flavour states a = e, /i, r to the eigentriad j at the scale ^/s. Explicitly, 
for the one-photon exchange diagram of Figure ^(a), we have: 



\pj Pj) 



(8) 



where s and s' denote the spins of the incoming and outgoing lepton i and 
r and r' those of the antileptons i^. 




(a) (b) 

Figure 3: Feynman diagram for the diagonal amplitudes. 

To actually evaluate (§), we have still to specify the values of the momenta 
Pj and p'j entering there. The point is that the amplitude for the two-body 
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reaction (|^ is of course a function of only two variables, which we may take 
to be the standard Mandelstam variables s and t, so that all components 
of the various momenta appearing in must be expressible in terms of 
them. The reasoning required to arrive at these expressions is not entirely 
trivial, but following the considerations given in which apply as well 
to the present case with but minor modifications, we obtain the following 
relationships between the different momenta to be used later for deriving the 
required expressions: 



Pj = ajq + bjq' + CjPi 
P'j = ajq' + bjq + Cjp[, 



(9) 



where 



{s — m'j + m?Y + stQ 
[s — ni?Y + ^t ' 
Cj{s + m^) — (s — m?j + m^) 

to ' 
Cj{s + + to) - (s - m| + + to) 
to ' 



(10) 



with to = t — 4m^, m the positron mass, and mj put equal to zero for reasons 
to be made apparent. 

In this paper, we shall be interested mainly in unpolarized cross sections 
which means that we shall need to evaluate sums of the absolute values 
squared of the amplitudes (|^) over all the spins s, r, s', r' . These spin-sums, 
as was explained in [0 , are not so readily performed as usual by the standard 
method of taking traces of 7-matrices because of the crossed terms between 
channels labelled by different j's obtained in squaring (|^). We shall therefore 
follow the tactics adopted in of explicitly performing the spin-sums in a 
specific Lorentz frame with a specific representation of the 7-matrices. As 
in [|r^, we choose to work in the cm frame of the channel i = 3 which in 
our convention denotes the mass eigenstate with the lowest mass m^, which, 
being in all cases considered at most of the order of the electron mass and 
therefore negligible, is put equal to zero. This gives: 



(E,0,0,u;), 

0, uj sin 6'^, oj cos 6'-^) 
(u;,0,0, -uj), 
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Pg* = (c<j, 0, — co'sin6'3, — ciJcos6'3), 
Pj = {Ej , 0, —Uj sin 9j , —uj cos 9j 
Pj" = {Ej, 0, —Uj sin 9j, —luj cos 6'^), (11) 



with 



2 



s — m 

u = 



2v^ 



5 



cos^^ = 1 + ^-^. (12) 



Further, from (^-(|rT|), one obtains 

Ej = ^^|2^[(s - m2)2 + st][{s - m] + rxi^f + sto] 

- (s + m2)[2(s-m2 -m^2) (13) 



with 



and 



a;, = jE|-m2, (14) 



sin 6',- = — ^ [is ^ vr?\\j is — vr?- ^ VTi^Y ^ stn 
2stQUJUJj I V ' ■' 

— (s — + m^) a/ (s — m^)^ + st 



sin 6'' = — — ^|(s + m^)^/(7^^mf~+m2)2~+st^ 
ZsiQUJUJj L V 

- (s-m^2 + m2 + to)V(s-m2)2 + st|, (15) 



with 



0^ = 0,+^;. (16) 

These formulae ([T2|)-([T5D, together with the formulae in the last paragraph, 
all reduce to the corresponding formulae derived in [0 for photo-transmuta- 
tion if we put the mass m of the e"*" in (^ equal to zero, which will indeed be a 
very good approximation in most applications. We have kept the dependence 
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on m explicit only for the sake of generality in case the formulae arc to be 
applied in future to other circumstances, such as lepton transmutations in 
/x^/x" collisions. 

We choose again for 7-matrices the Pauh-Dirac representation: 
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1 
-1 



; 7 



O-fc 
-CTk 



(17) 



The spins of the incoming and outgoing leptons j we quantize along the di- 
rection p3 and P3 respectively, while the spins of the e"*", whether incoming 
or outgoing, we quantize along its direction of motion. With these specifica- 
tions, the wave functions of the leptons j are given by: 



u+{pj) = 



^2{Ej + mj) 



( 



I + rrij \ 








'j 








(18) 



and 



2^2{Ej + mj) 



2^2{Ej + mj) 



( (i?,+m,)(l-e-'^3) X 
(£;,■ + m,-)(l + e-^3) 

( (E^ +mj)(l + e 
(Ej- +mj)(l - e 



V 



) 



(19) 



while for the e"*", we have the wave functions: 



/ \ 



-ijj 



m] 



\ E + m J 
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v-iq) 



^2{E + m) 



UJ 



E + m 




(20) 



and: 



v4q') 



2^/2{E + m) 



-ie'- 



( -07(1-6-^^3) \ 

-cu(l + e-'^3) 
(E + m)(l - 

V (E + m)(l + e-^^3) ) 

( cu(l + e-^^3) \ 
cu(l - e-^^3) 
(E + m)(l + e-*^3) 

V (E + m)(l y 



(21) 



With the wave functions in (p^)-(|2T|), it is straightforward to evaluate 
the diagonal amplitudes in (^. We obtain the following: 



{M 



ie"^ r E 

- — I j(l + cose^)[(E, + m,) + (E^ - m,)e- 



+co'ti;j(l — cos6'3)e 



-id. 



(1 + COS 613)6 



-id. 



%e 



:^(l + cos^^){|[(i?,+m,) + (E, 



UJUJi 



2 

3 
2 



^^(1 - cos^^){(E, + m,) - (E, - m,)e-2^^^ } 



3 > — 



-- sin^^^^ (E, + m,) + (E, - m,)e-^^^^ 



-^|sin^^{(E, +m,) - (E, -m,)e-2*^^}, 



(22) 
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where subscripts denote the spins of the leptons j and superscripts the spins 
of the e"*", with the right index pertaining to the incoming and the left index to 
the outgoing particle. Using the formulae derived earlier in (|T2|)-(|T3p, these 
amplitudes can then be expressed in terms of the Mandelstam invariants s 
and t as desired. 

3 The Reaction Amplitude (b) 

Turning next to the diagram of Figure [l|(b) which contributes to all three 
reactions we proceed in a similar manner. The transmutational 

amplitude for the reaction: 




(23) 



is given by a rotation in generation space, thus: 




(24) 



from the diagonal amplitudes Aij depicted in Figure ^(b) for the reaction: 

e+e- ^ ipj (25) 



where 



{Mf^Ys': = ^e^[^.'(p;)7'^^..'(gp]^— ^[^.(g)7M«.(p)]- (26) 




We work now in the cm of the incoming and e system with 




and 



E = Ej = v^/2, 



cos ^3 




(28) 




2ti;co'o 
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Further, with the spins of the outgoing particles quantized along qi^ 

and P3 respectively, and those of the incoming e^,e~ along their directions 
of motion q and p respectively, the wave functions are given by: 



u+ip) 



u-{p) 



v-{q) = 



^J2{E + m) 



^/2{E + m) 



m) 



( \ 

E ^ m 


( E + m\ 


—UJ 




—UJ 



\ E + m J 

CO 



E + m 




and 



2^2{E,+m,) 



2^2{Ej + mj) 



2^2{Ej + mj) 



( (Ej +mj)(l - e 

/ (£;,+m,)(l + e-^^ 
[E,+m,){l-e--"'^ 
— + e 

(Ej+mj)(l-e 
V (^, +m,)(l + e 



-ie'y 



J 



(29) 



) / 
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Hence, one obtains the amphtudes: 



( a;j(l + e 

(E,- + m,-)(l + e-''^ 
V (E,+m,)(l-e-^^^) / 





: = (Mf>)z 


= mm,- cos ^w, 






+ = EEJl + cos9'^ 

s 




: = (<): 


+ a' 
T = — mrrij cosu,, 

~^ s 


{Mf)t 


; = Off): 


^ = — Errij sin 6'', 
s 






_i_ = Enij sin 6'^ 


{MfY- 


: = Off); 


i = — mE",- sin 




t = (A<<")1 


; = EEj{i-co^d'- 


{MfY-- 


: = (-Mf); 


= mEj sin 6'', 



(30) 



which again can a 
and t by means of the formulae in 



(31) 



1 be expressed in terms of the Mandelstam invariants s 



4 Spin-summed Differential Cross Sections 

Substituting the diagonal amplitudes in ([2^ ) and (|3TD into respectively the 



formulae (y) and (|2^) and adding the two contributions, one obtains the spin- 
amplitudes for the actual transmutation reaction (H), (|) and (§). Hence, 
taking the absolute values squared of these amplitudes and summing over all 
spins, one obtains the spin-summed differential cross sections as desired: 
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for energies measured in GeV, where uj and uo' are respectively the cm mo- 
menta of the actual incoming and outgoing particles in the transmutation 
reaction. Although this is in principle straightforward, a few practical obser- 
vations are in order. 

First, in performing the sum in and (^) over the diagonal states j, it 
is convenient to make use of the unitary property of the rotation matrix Saj 
to write, for a ^ 13: 



S^,SljM, = S^iSli\M, - M^l (33) 



which, as explained in |T^, is a good approximation whenever the mass 
eigenvalues rrij are hierarchical and avoids the need to know the rotation 
matrix to unreasonably high accuracy. Besides, as we shall see, it gives us 
a clearer picture of how transmutation cross sections behave as functions of 
the Mandelstam invariants s and t. 

Second, as can be seen in (|33|), the transmutation amplitude being pro- 
portional to A^i — with the two amplitudes differing just by the mass 
values, i.e. whether mi or ms, the cross section for transmutation is at most 
of order m\/ s and decreases rapidly with increasing energy. For high s, there- 
fore, M.I and M.^ will largely cancel leading potentially to inaccuracy in a 
direct calculation with the formula (|33|) . Indeed, this was exactly what we 
found in our actual calculations, especially in the the amplitude (a) where 
other large cancellations occur in the exact formulae. This computational 
difficulty can be avoided just by expanding the amplitudes to order m\/s 
giving: 



s + t 



5 



ie^ , , rm? .m? / —t 



t I 2s s \ s + t 



js + t (nil .iTT'i 



2 y —t \ ^/s s y s + t 



(34) 



and all other components zero, where we have also neglected terms of the 
order of the electron mass m. This approximation is already very good by ^/s 
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of order 10 GeV, at least near the forward direction where the amphtudes are 
large, and becomes eventually necessary above this energy for computations 
without double precision. In Figure ^ is shown the spin-summed differential 
cross sections for the reaction e~^e~ — > e~^T~ at ^/s = 10 GeV, 100 GeV 
calculated with the rotation matrix element Sai of the DSM scheme taken 
from ref. [|l^, The curve at 10 GeV is calculated with the exact formulae 
dH) which is seen to be almost indistinguishable from the crosses calculated 
with the approximate formulae (0). The curve at 100 GeV is calculated 
with (^4|) where the exact formulae is found to have problems with accuracy 
in application. 

Third, we note that the two sets of diagonal amplitudes (^) and (|3lD 
were each calculated in a particular Lorentz frame, namely for the diagram 
(a) in the cm frame of the e"'"£j system and for the diagram (b) in the 
cm frame of the incoming e^e~ system. Although the amplitudes were all 
converted in the end into functions of the invariants s and t, the directions 
of spin quantization are still frame- dependent. Hence, strictly speaking, the 
two frames for (a) and (b) being different, the two respective spin- amphtudes 
could not be added in the manner that we have done above. However, the 
electron mass m is so small compared to the energies we are interested in that 
this difference in frame is entirely negligible for practical purposes. Were the 
present formalism to be adapted in future to say fi~^fi~ collisions, then this 
would be a point to be borne in mind. 

With these points clarified, we have not encountered any more practical 
difficulties in computing the cross sections of the three transmutation reac- 
tions (|D-(@). Rather than presenting our results in a wide range of s and t, 
which could be confusing, we shall instead first give here a description of the 
general features, and then in the next section a detailed report on the result 
at = 10.58 GeV, namely at the T(4^) where BaBar [P, Belle fH, and 



Cleo 1 23] have already collected a massive amount of data, in principle ready 
to be confronted with the above predictions. 

Consider first the reactions (H) and which receive contributions from 
both diagrams (a) and (b), and are very similar except for the difference 
in normalization due to the different values of the rotation matrix elements 
Sal- As in ordinary Bhabha scattering, the amplitude (a) is dominated by 
the pole at t = which gives the cross section a sharp forward peak, as can 
be seen in the examples of Figure ^ Except at large scattering angles where 
t is of order s, this peak overshadows the contribution from the (b) diagram. 
However, the forward peak for the transmutation reactions and is 
nowhere near as sharp as for ordinary Bhabha scattering, as can be seen in 
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Figure 5: The ratio at 10 GeV of the cross section of reaction e~^e — > e^r 
over that of ordinary Bhabha scattering e+e" — > e+e" as a function of t. 



Figure ^ The reason for this difference is seen in (|3^), where one notices 
that the normally dominant spin non-flip amplitudes with 1/t behaviour are 
of order m\/s, while the spin flip with a weaker behaviour are of 

order m\j The same formulae p^ ) explains also the sharp decline of the 
cross section with increasing energy as well as its change in t-dependence as 
the spin flip terms become ever more dominant, both of which effects can be 
seen in Figure H by comparing the curves at \fs = 10 and 100 GeV. 

The other reaction (H) receives contributions only from the (b) diagram 
which has no peak in the forward direction. It is distinguished from the 
same diagram in ordinary Bhabha scattering by the fact that, like the (a) 
transmutation amplitude, it is also dominated by the spin flip terms at high 
energy. Without the sharp singular peak in the forward direction, it gives, 
in contrast to reactions (@) and @, a finite total cross section, the rough 
energy dependence of which is shown in Figure where one sees that, as in 



photo-transmutation [17|, the cross section rises shortly after threshold to a 
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Figure 6: Cross section for the reaction e+e 
scattering angles as a function of y/s. 



fi~^T integrated over all 



peak and then declines as -y/i increases. 

We note that since the rotation matrix elements Saj enter only in the 
normalization of the cross section and are themselves only weakly dependent 
on energy, their actual values do not affect the discussion above on the the t- 
dependence, and qualitatively also on the s-dependence, of the cross section. 



5 Transmutation Cross Sections at ^/s = 10.58 
GeV 

The reason for selecting this particular energy corresponding to the mass of 
the T{4S) for detailed analysis is that two experiments of ultra-high sensi- 
tivity, namely BaBar and Belle, are running and have collected already up to 
20 fb~^ of luminosity, with much more expected in the near future pT| , 
Although these experiments were designed originally to look for other rare 
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effects like CP-violation from B decay, their data could conveniently be used 
also to search for the transmutation reactions of interest to us here. 

Consider first the reaction (|^) with contributions only from Figure |I](b). 
The normalization of the cross section is, according to (^), given by the 
rotation matrix elements 5*^1 5*^1. If one accepts the contention of that 
it is the rotation which is giving rise to fermion mixing and mass hierarchy, 
then these rotation elements should be related to the data on the mixing pa- 
rameters and mass ratios. Indeed, as explained in that paper, the quantity 
SfiiSri would then be given to a good approximation by sin 9 cos 9 where 9 is 
the rotation angle from the r-mass scale where the mass matrix is diagonal in 
the lepton flavour states to the scale 10.58 GeV of say the BaBar experiment. 
Its value can thus be read off from Figure 3 of |]15[, either directly by inter- 
polating the actual data points or equivalently, since the rotation curve from 
the DSM scheme is seen there to be a very good fit to the data, by taking 
the values from the DSM calculation, in either case leading to an estimate 
Sfj.iSri ~ 0.043. With this then for the normalization, one obtains Figure 
for the spin-summed diffrential cross section of reaction (^). Integrating over 
the whole angular range, one obtains a cross section of around 80 fb. This 
can also be read in Figure ^ which was in fact calculated already with this 
normalization. In practical terms, this could mean as many as 1600 events in 
the data sample of 20 fb~^ already collected by BaBar, assuming 100 percent 
efficiency. 

Consider next the reaction which receives contribution from both 
diagrams in Figure |I|. Its normalization depends now instead on the rotation 
matrix elements SeiSri which is considerably smaller and harder to estimate 
directly from data. Let us then just insert the values obtained from the DSM 
scheme pO| which was shown to fit the data very well and can be taken 
as a convenient interpolation device. The actual value can be read in Figure 



4 of at 10.58 GeV to be around 0.0092. This then gives Figure g for 
the spin-summed differential cross section for this reaction. As in ordinary 
Bhabha scattering, the cross section for (H) is divergent at t = 0. However, 
this region cannot be explored experimentally, the detectors being insensitive 
to the forward region with |t| < 5 GeV^ A rough estimate from Figure 
^ then yields for the integrated cross section for reaction (Q) over the range 
\t\ < 5 GeV^ a value of about 20 fb, which means about 400 events in the 
sample of 20 fb~^ already collected in last year's run assuming again 100 
percent detection efficiency. Compared to reaction (^) above, we note the 
very different angular dependence which is here dominated by the diagram 
(a) of Figure ^ absent in reaction (^). The reason why the cross section for 
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is smaller than that for despite the fact that (|^) receives contribution only 
from the sub-dominant diagram (b) is the very different sizes of the rotation 
matrix elements involved, (5'^i/5'ei)^ being ~ 20, as read from Figure 4 in 

However, we should perhaps stress once more that the Figures ^ and ^ 
calculated by the method of fl^ represent only the kinematic effects of the 
rotating mass matrix. In particular, though calculated with DSM rotation 
matrices, the transmutation cross sections given are not those predicted by 
the DSM scheme in which, as seen in |jl9l, the above kinematic effects are 
largely cancelled by other rotation effects implied concurrently by the renor- 
malization mechanism driving the mass matrix rotation. 

Apart from minor kinematic differences, the spin-summed differential 
cross section for the reaction @ has the same t-dependence as (H) but a 
different normalization, namely with Sri in (0) replaced by S^i in (^. For 
instance, taking again the rotation matrix elements from Figure 4 of |1T7| , 
one obtains a cross section for (^) a factor ~ 2 x 10~^ smaller than for (^, 
making it probably difficult in any case to observe in the near future. 

For all 3 reactions, there is in principle also a contribution from the trans- 
mutational decay of the T(4S') resonance with this mass, but this will be seen 
in the following section to be negligible in comparison with the above contri- 
butions. 

From the above results, it would appear that if one accepts the interpreta- 
tion as given in |T^ of fermion mixing and mass hierarchy as rotation effects, 
and that there are no other rotation effects than that of the mass matrix, then 
there will be lepton flavour-violating transmutation effects in e'^e~ collisions 
at BaBar, Belle, and Cleo energy, which are of a magnitude to be observable 
by these experiments. On the other hand, if the interpretation of JlSl is not 
accepted, one knows at present of no other empirical means for estimating 
the rotation angles and thus no estimate for the absolute rates of transmu- 
tation can yet be made. However, since the rotation matrix elements enter 
only in the normalization of the cross section, the calculation above is easily 
adaptable to any rotation matrix obtained from whatever source, whether 
empirical or theoretical. 

Suppose that the reactions and are indeed observed, can one be 
sure that they are due to transmutation and not some other lepton- violating 
effect? The answer to this question would seem to be quite affirmative since 
one has here the differential cross sections as functions of 2 variables, each 
with distinctive characteristics. For example, it is predicted that the cross 
section for (^) should be peaked sharply forwards as seen in Figure H, while for 
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(|^) it should have a roughly sin^ 9 behaviour, as seen in Figure ^ And both 
are predicted to be spin-flip dominated, which assertion may be verifiable 
with the decaying r serving as its own spin-analyser and the spin- dependent 
cross section calculable from the amplitudes (P^ ) and when the occa- 
sion demands. Thus, if the reactions are observed at all with any reasonable 
statistics, the signatures for transmutation would seem to be quite unmis- 
takable. 

We have restricted the discussion in this section to only the operation 
energy of BaBar, Belle and Cleo, but very similar remarks apply also to the 
BEPC energy ||2^. The expected transmutation cross sections at BEPC, as 
seen in Figure ^ for reaction (^, are even larger, but this advantage is un- 
fortunately more than offset by the lower luminosity so far achieved by the 
machine. For this reason, for BEPC, we think a search for transmutation 
in if) decay will be more immediately profitable, as we shall elucidate in the 
following section. As for LEP, according to Figure ^ for example, transmuta- 
tion cross sections would have fallen much below the fb level by that energy 
even for the optimistic scenario of [|1^] and are thus unlikely to show up in 
the data collected. 



6 Vector Boson Decay 

Vector boson formation and decay occur naturally in e^e~ collision so that 
they will have in any case to be taken into account in studying transmutation 
from this process. Besides, being essentially a single particle effect, vector 
boson decays are easier to analyse theoretically and can give more succint 
conclusions than the reactions previously considered. We propose therefore 
to examine in this section the following transmutation decays of vector bosons 
into leptonic final states: 

V ijp, (35) 
for £a(£/3) being e,/x or r. Order of magnitude estimates for some of these 



decays have already been made in |TB[. Now, with the procedure developed 



in |]T^ an explicit calculation can be carried out. 

Consider then in general a decay of the type (|35D as depicted in Figure 
11(a). The decay amplitude is to be evaluated at the scale y/s = M, with M 
being the mass of the decaying boson. At this scale, the fermion mass matrix 
m is in general not diagonal in the fiavour states e, fi, t but in some other 
states, say, i = 1,2,3 with masses (eigenvalues of m) m^, the two triads of 
state vectors being related by a rotation matrix Sai = {a\i). To calculate 
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the decay amplitude, we first evaluate the amplitudes for the decays into the 
diagonal states i, namely V — > £iii as depicted in Figure 0(b). For a vector 
boson with polarization vector e^, this amplitude is: 

Mi = ge^Sk)ur{pi)-if,Vs{qi), (36) 

where r and s denote the spin respectively of and Ii. We do not need 
to specify the coupling strength g of the vector boson to the lepton pair 
for in our branching ratios calculation it will be cancelled out. According 
to the procedure suggested in the amplitude for the decay of interest 
([35| ) is then given just by a rotation of these diagonal amplitudes M.^ to the 
appropriate lepton flavour states, thus: 

M"-^ = Y.SlM.Sp,. (37) 

i 

To obtain the decay rate, with which alone we shall be concerned at 
present, we average over the initial polarization a of the decaying boson and 
sum over the final spins r, s of the product leptons, obtaining: 

a,r,s r,s r,s 

where 

-A^^fs = 9^SliUr{Pihfj,Vs{qi)Sf3i, (39) 
i 

-^S- = 9J2^^^^''rip^)Hiq^)Sp^. (40) 

i 

The remaining sum over lepton spins again cannot readily be done as an in- 
variant trace because of the crossed terms between different internal channels 
i. We proceed as for the reactions (0) — (^ above. 

We choose to work in the rest frame of the decaying vector boson: 

k = + = 0, (41) 

so that: 

pf = (E,0,0,-^,), 

gf = (^,0,0,^,), 

fc^ = (M, 0,0,0), (42) 
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with 



E = M/2, 



jt;2 _ 
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As a result, we have 
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With these exphcit expressions for the lepton wave functions, it is easy to 
evaluate the spin amplitudes and perform the sum over the lepton spins 
giving: 

3 ^ 3 



:/3 |2 
ars I 



5*^4 m j m j S*^ ^- 5*0, j . 



(45) 



The states i being by definition the eigenstates of the mass matrix m at the 
scale of the decaying boson mass M, and Sai the rotation matrix relating 
these states i to the lepton flavour state a = e, /i, or r, the sum S^^iSfuvtii 
is just the element {pi\m\(5) of the matrix m at the boson mass scale. Hence, 



q / J \-' 'cm 



,fi |2 
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This gives the total width for the decay V 



i(i as: 
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where we note that the phase space factor ul^p/ E^Ep refers as per ||T3 to the 
freely propagating "external" leptons, with momentum and energy respec- 
tively: 
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More conveniently, since the coupling g has not been specified, one writes for 
a P, i.e. transmutational decays: 



2 



\{a\mm', (49) 



where we have neglected terms of the order of the electron mass me compared 
to M. Multiplying then this ratio by the experimental branching ratio, if 
known, of the boson V decaying into e^e~ gives the branching ratio of the 
transmutational £a^/3 mode. 

One notices that apart from a numerical factor of 2 and the phase space 
factor uj"^^/ EaE/s, the formula for the branching ratio (^) for transmuta- 
tional decays is the same as the order-of-magnitude estimate given in |TB 



This formula is supposed to hold in general and can be applied to calcu- 
late the branching ratio of transmutational decays of the type ( |55D given the 
matrix element (a|m|/3) at scale M. 

For a numerical example, we take again the interpretation in of exper- 
imental data to estimate the rotation matrix elements required or, in case the 



data is insufficient, employ the DSM result of |2^ as an interpolation formula, 
for which the required mass matrix elements {a\m\l3) can be read directly 
from Figure 3 of [|16|. With these elements, the formula (^) can be applied 
immediately to calculate the branching ratios of transmutational modes in 
the decay of any vector boson by normalizing to the empirical branching ra- 



tios of the e^e~ mode (or to the mode if more accurate) given in [|T4 

The result for the most experimentally interesting vector bosons listed in 
is given in Table ||. (Again, the reader is reminded that these numbers 



represent only the transmutation effects due to kinematics of the rotating 
mass matrix alone but not, though calculated with DSM rotation matrices, 
predictions of the DSM scheme in which cancellations with other rotational 
effects occur [|19|.) We have calculated also the branching ratios for the other 
vector bosons listed in WM, such as the higher excitations of ip, but for these. 



once the mass gets above the DD threshold, hadronic decays prevail, leading 
to large total widths and hence uninterestingly small branching ratios for the 
modes that concern us here, as can be seen in the example given of ?/'(3770). 
In any case, the predicted branching ratios for all higher excitations of and 
T are very similar to those given in Table [l| for ?/'(3770) and T(10860) and 
are therefore not given again there. 

In Table |l], we note first that the predicted branching ratios for fiavour- 
violating transmutational decays even for the fast rotation implied by [|15| are 
not as large as one might fear at first sight. Indeed, all the estimates survive 
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Boson 


Mode 


Predicted Branching Ratio 
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ZD 1 1 
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1 1 
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1.9 X 10-^ 




er 
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fJ>T 
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er 


8.8 X 10-9 
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3.5 X IQ-^-J- 



Table 1: Predicted branching ratios for transmutational decays. Modes which 
could be accessible to experiment in the near future are underlined. 
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existing experimental bounds comfortably, which are surprisingly weak for 
vector boson decays. Of the vector bosons listed, the data book |]T1| gives 
upper limits on lepton-fiavour violations only for Zq decay, which for the 
modes fiT,eT, and e/i are respectively 1.2 x 10"^, 9.8 x 10~^, and 1.7 x 10~^ 
in branching ratios, which are seen to be all easily satisfied by the estimates 
in Table ^ 

Secondly, we note several entries of branching ratios in Table |l] fall well 
within the sensitivity range of present experimental set-ups. In particular, 
the fiT decay modes of '0(15'), '0(25'), T(IS'), and T(25), are seen each to 
have a predicted branching ratio of several parts in a million and hence well 
within the present sensitivity range of BEPC [Q, PEP II (BaBar) and 



BELLE ||2^. For instance, BEPC has already collected more than 20 million 
■^'s to-date, and expects to collect twice as many more next year, which would 
mean for the predicted branching ratio of 6.3 x 10^^ as many as 120 now, 
and 360 next year, of fir decays assuming 100 percent detection efficiency. 
Similarly, PEP II has already accumulated for BaBar in a year over 20 fb~^ 
in luminosity at T(45), which means that if the machine is run at the T(15) 
with a cross section of around 25 nb, it would collect in just a couple of 
months of running already enough T(15) events to give, at a branching ratio 
of 1.2 X 10~^, over 100 fj,T decays, again assuming 100 percent detection 
efficiency. We note in passing, however, that at the energy 10.58 GeV at 
which the present BaBar experiment is run, the T(45) resonance at that 
mass, being above the BB hadron threshold, has a branching ratio of only 
2.8 X 10~^ into the already most copious transmutational fir mode, which 
makes the resonance effect negligible compared with direct transmutation via 
the process of Figure |l|(b) as calculated in the preceding section. 



7 Concluding Remarks 

In summary, our conclusions are as follows. Although flavour conservation 
has been checked to very high accuracy in leptonic decays such as — *• ej 
and /i — > eee, it may be premature to conclude that flavour-violation will 
always be small. The fact that the fermion mass matrix can rotate with 
changing scale can mean that flavour though accurately conserved at some 
energies are appreciably violated at other energies. We suggest therefore that 
lepton flavour violation should be routinely tested by experiment whenever 
conditions are favourable. Any sizeable signal (e.g. BR in vector boson of 
say order greater than 10^^^ ^5|) can mean mass matrix rotation, and any 
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information on rotation, and especially on its energy dependence, can give 
us new insight into the problem of fermion generations. 

If we make the assumption that fermion mixing and mass hierarchy are 
due to mass matrix rotation, for which hypothesis there seems to be some 
empirical support [|T^], then the magnitude of the kinematic effects on flavour- 
violating from a rotating mass matrix can be estimated, and have been found 
to be appreciable and very likely detectable with the present generation of 
high sensitivity experiments such as BaBar, Belle, Bepc and Cleo. A detec- 
tion of the effect at the estimated level would lend support to the hypothesis 
of |jl5| which would be a big step forward towards the solution of the gen- 
eration puzzle. However, a negative result, unfortunately, cannot rule out 
the rotation hypothesis, for it can happen, as it is shown to be actually the 
case in a specific example in |T^, that the effect of the rotating mass ma- 



trix calculated here is modified or even cancelled by other rotation effects 
accompanying the mass matrix rotation. 

We thank Gian Gopal and John Guy for useful discussions on the exper- 
imental feasibility of testing the ideas of this paper and for their encourage- 
ment through all the ups and downs of its progress. 
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